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Abstract. A rationality result previously proved for Robertson-Walker met¬ 
rics is extended to a homogeneous anisotropic cosmological model, namely the 
Bianchi type-IX minisuperspace. It is shown that the Seeley-de Witt coeffi¬ 
cients appearing in the expansion of the spectral action for the Bianchi type-IX 
geometry are expressed in terms of polynomials with rational coefficients in 
the cosmic evolution factors wi(t), W2(t), ws(t ), and their higher derivates with 
respect to time. We begin with the computation of the Dirac operator of this 
geometry and calculate the coefficients ao , a2 , <24 of the spectral action by using 
heat kernel methods and parametric pseudodifferential calculus. An efficient 
method is devised for computing the Seeley-de Witt coefficients of a geometry 
by making use of Wodzicki’s noncommutative residue, and it is confirmed that 
the method checks out for the cosmological model studied in this article. The 
advantages of the new method are discussed, which combined with symmetries 
of the Bianchi type-IX metric, yield an elegant proof of the rationality result. 
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1 . Introduction 


Quantum cosmology studies the early universe where the energy scale is so high 
that one would need to incorporate into the theory both quantum gravity and the 
Standard Model with unbroken symmetries such as the hypothetical supersymme¬ 
try. This makes the exact solution of any quantum cosmological problem essentially 
impossible as it amounts to solving quantum fields genuinely interacting with quan¬ 
tized gravity, the nature of which we know little about. However, since we are only 
interested in the large-scale behavior of the early universe whose mass-energy dis¬ 
tribution is highly homogeneous, we are encouraged to exploit this symmetry and 
focus only on a few long-wavelength degrees of freedom. In quantum cosmology, 
this common practice is known as the minisuperspace approximation, which can be 
rigorously justified under certain criteria [521 35] , 

Well-known examples of minisuperspaces include the Robertson-Walker model 
and its anisotropic generalization to the Bianchi type-IX model. The Robertson- 
Walker metric is of the form 


ds — dt~ a^t') da~, 


where a(t) is a general cosmic factor of the expanding universe and da 2 is the round 
metric on the 3-sphere S 3 . The Bianchi type-IX model, which enjoys a reduced 
SU (2) isometry group rather than the full § 3 symmetry, is written as 


(1) ds 2 = VJi(t)w 2 (t)w 3 (t)dt 2 


W 2 (t)w 3 (t) 2 w 3 (t)w 1 (t) 2 W 1 (t)w 2 (t) 2 

-— CT 1 - — a 2 H-7T\— 03 > 

Wl{t) w 2 (t) w 3 (t) 


where eq are left-invariant 1-forms on SU( 2)-orbits. 

Beside trying to fit the Standard model into the early universe spacetime, another 
possible investigation of physical importance is to quantize supersymmetric systems 
in a Bianchi-type IX minisuperspace. Since the energy scale of the early universe 
is very likely to be higher than the supersymmetry breaking scale, it is of great 
interest to see how supersymmetry may possibly change the picture of the early 
universe. Since the Hamiltonians of supersymmetric systems can often be identified 
with well-behaving elliptic operators acting on spin bundles for Majorana fermions 
or fields of differential forms for Dirac fermions, the spectral action [5] approach 
appears to be a handy tool, where the celebrated Atiyah-Singer index theorem may 
be used to calculate the Witten indices of certain supersymmetric theories and thus 
to conclude whether the supersymmetry can be spontaneously broken 02112]. 

Nonconmrutative geometry m describes geometric spaces by spectral triples 
(A, "H,!?), where A is an involutive algebra represented by bounded operators on 
a Hilbert space 7i , and D is an unbounded self-adjoint operator in T-L that plays 
the role of the Dirac operator by encoding the metric information. This set up, 
which includes a great variety of nonconnnutative spaces, generalizes Riemannian 
geometry since Connes’ reconstruction theorem m states that if A is commuta¬ 
tive, then, under suitable regularity conditions, the triple consists of the algebra of 
smooth function on a spin c manifold acting on the L 2 -spinors and D is the Dirac 
operator. 

Given a spectral triple, the spectral action principle mm considers the func¬ 
tional, 

Trace(/(D/A)) ~ £ / \D\~? + /(O)Cd(O) + • • • , 

/3e n ^ 
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as the fundamental action functional, where / is a positive even function on the 
real line and A is a real parameter that fixes the mass scale. The details of the 
above asymptotic expansion are explained in m- For a spin c manifold M, the 
coefficients of this expansion are determined by the Seeley-de Witt coefficients a n 
appearing in the heat expansion, 

OO 

Trace(e -tI?2 ) ~ *- dim ( M )/ 2 ^ a 2 n (D 2 )t n (t -)■ 0+). 

n=0 

In addition to extracting local geometric information and recovering the Einstein- 
Hilbert action, the spectral action generalizes the theory to include noncommutative 
spaces which give rise to the Standard Model gauge fields with Yang-Mills actions as 
well as couplings to the modified gravity [51131 IT31I551 HD] . In this sense, this principle 
is able to unify the Standard Model with general relativity under a generalized 
notion of Riemannian geometry. 

Considering the variety of noncommutative geometric spaces and the physical 
implications of the spectral action MEHU], it is of great importance to use and 
develop different methods for computing this action. In particular, any interpreta¬ 
tion of the full expansion of a spectral action is desirable. For example, using the 
Poisson summation formula, the spectral action for the Dirac operator on highly 
symmetric manifolds, such as products of spheres by tori, was computed in [8|. 

A considerable amount of work has been carried out on the cosmological impli¬ 
cations of the spectral action in recent years [32 IS3 [Ml [Ml @01112 [ 12 031 [IB] ■ 
For the Euclidean Robertson-Walker spacetime with a general cosmic factor a(t), 
Cliamseddine and Connes have devised an efficient method in [9] for computing 
the terms of the spectral action, which is based on making the use of the Euler- 
Maclaurin formula and the Feynman-Kac formula. They computed the terms up 
to aio in the expansion and made a conjecture, which was addressed in [21] by us¬ 
ing pseudodifferential operators and heat kernel techniques. That is, it was shown 
that a general term in the expansion is described by a polynomial with rational 
coefficients in a(t) and its derivatives of a certain order. 

The present paper is intended to study the spectral action for the Bianchi type- 
IX minisuperspace, which is a homogeneous anisotropic cosmological model. Since 
the spectral action of a geometry depends on the eigenvalues of the square of its 
Dirac operator, we explicitly compute the Dirac operator D of the Bianchi type-IX 
metric in f|2] and derive the pseudodifferential symbol of D 2 . In fj3l following a 
brief review of the heat kernel method that uses pseudodifferential calculus for the 
computation of the Seeley-de Witt coefficients [2F,, we present the calculation of 
the terms 00 , 02,04 in the expansion of the spectral action associated with D. 

We devise a new method for calculating the Seeley de-Witt coefficients of a 
geometry in terms of noncommutative residues of operators, which extends the 
result on the realization of the Einstein-Hilbert action as the residue of a power 
of the Laplacian, see [281 SOI E5]. This method is explained in detail in i[4] and 
it checks out to give the same result for the calculated terms 00 , 02,04 for the 
Bianchi type-IX metric. Combining the symmetries of the metric with technical 
properties of pseudodifferential symbols of parametrices of the Laplacians, which 
significantly simplify in view of the new method using the Wodzicki residue, we 
prove a rationality result for a general term in the expansion of the spectral action 
for the Bianchi type-IX metric in S[5] That is, we show that general terms of the 
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expansion are expressed by several variable polynomials with rational coefficients 
evaluated on W\ (t ). W 2 (t). W 3 (t ), and their derivatives of certain orders. In 33 we 
discuss the gravitational instantons, which form an especially interesting class of 
Bianchi type-IX models, and elaborate on the significance of the rationality for 
the spectral action in relation to the arithmetic and number theoretic structures in 
mathematical physics. Finally, our main results and conclusions are summarized 

in 33 


2. The Dirac Operator of Bianchi Type-IX Metrics 

The heat kernel method that uses pseudodifferential calculus for computing the 
Seeley-de Witt coefficients of an elliptic operator on a compact manifold relies on the 
pseudodifferential symbol of the operator in local charts )26j . Thus, in this section 
we compute the Dirac operator D of the Bianchi type-IX metric and thereby obtain 
the symbol of D 2 . 

The most efficient way of computing the Dirac operator of a geometric space is 
to use an orthonormal coframe {9 a } for the metric, which, from the definition of 
D, yields 

a 

where V s is the spin connection of the spin bundle and {9 a } is the predual of 
the coframe, see [24]. Since V s is the lift of the Levi-Civita connection V to 
the spin bundle, one starts with computing the matrix of 1 -forms w = (cv b ) such 
that V = d + oj in terms of 9 a in a local chart (x^), which can be lifted to the 
matrix of the spin connection 1-forms by making use of the Lie algebra isomorphism 
/r : so(m) —> spin(m) given by 

M = - a Y, A ab e a e b , A = ( A ab ) e so(m), 

a,b 

where m is the dimension of the manifold. We note that {e a } is the standard 
basis for R m considered inside the Clifford algebra of R m where spin(m) is spanned 
linearly by {e a e b ;a < b}. 

The u> b are found uniquely by writing 

V<9“ = 

b 

and imposing the conditions that characterize the Levi-Civita connection, namely 
metric-compatibility and torsion-freeness which respectively imply that 

<4 ? = -<*& d9 a = Y J u a b /\9 b . 

b 

Therefore the Dirac operator is written as 

d = X + \ X Y<-rV. 

a,b,c 

where uj b c are defined by 

c 

and the matrices 7 “ represent the Clifford action of 9 a on the spin bundle, namely 
that they satisfy the relations ( j a ) 2 = —I and 7 ° 7 b + 7 b 7 ° = 0 for a^b. 
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We go through this process for the Bianchi type-IX metric, 


j 2 ^ u , wiW 2 cos(r?) , wiw 2 cos(r?) 

as = W1W2W3 at at ^- dtp dip H- dip dtp 


w 3 w 3 

2/„\ -o2/ 

(Cl 


fu>2W 3 sin 2 (77) cos 2 {ip) | ^ 103 sin 2 (?j) sin 2 (7/;) i W2COS 2 (i]) 


\ Wi \ w 2 

(w\ — u>|) W3 sin(77) sin(i/0 cos(i />) 


w 3 


wiw 2 

(w\ — u>|) W3 sin(?7) sin(i/’) cos(?/>) 
UI1W2 


-dr/ dp 


dp drj 


(w 2 w 3 sin 2 (ip) u>iu; 3 cos 2 (^)\ W1W2 , ( ,, 

- 1 - drjdrj-l - dip dip, 

\ 


U> 2 


W 3 


d(j> d(p 


which is written in the local coordinates (x^) = (t, 77, <j>, ip), where S 3 is parametrized 
by the map 


(■ r],(p,ip) i->- ^cos(r// 2 )e*^ + ^^ 2 ,sin( 7 y/ 2 )e*^ ^)/ 2 ^ , 


with the parameter ranges 0 < 77 < 7r, 0 < p < 2 n, 0 < ip < 47r. An orthonormal 
coframe for ds 2 is given by 


0° 

9 1 

e 2 

e 3 


y/W\W2W 3 dt, 

sin(77) cos(^) 
sin(77) sin(-0) 
cos(??) 


w 2 w 3 
■w 1 

W1W3 

w 2 


dcp — sin(i/>) 


dtp + cos(ip) 


w 2 w 3 

Wi 

WiW 3 

U>2 


dr], 


dr], 


W1W2 


w 3 


dtp - 


wiu;2 

w 3 


dip. 


By explicit calculations in this basis we find that the non-vanishing co b c are 
determined by the following terms: 


Wn = -- 


w 33 — “ 


w 32 _ 


VJ2 {WiW 3 — w 3 w[) + WiW 3 w' 2 
2 (w\W 2 W 3 ) 3 / 2 
W2 (w 3 w[ — W\W 3 ) + W\W 3 W 2 
2 (w\W 2 W 3 ) 3 / 2 
w\ (w| - wf) + wlwl 
2 {w\W 2 W 3 ) 3 / 2 


w 22 — — - 


w 23 — 




W2 (w 3 w[ + W\W 3 ) — W\W 3 W 2 
2 (w\W 2 W 3 ) 3 / 2 
w\wl - wl (wl + wf) 
2 (wiW 2 U> 3 ) 3 / 2 
w 2 wj - w\ {wl + wf) 
2 {w\W 2 W 3 ) 3 / 2 
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Thus we achieve an explicit calculation of the pseudodifferential symbol of the Dirac 
operator which is written as 


a(D){x,0 = ^* 7 a e^ + i + 



rfy/wi (csc( 77 ) cosQ) (g 4 cos(77) - g 3 ) + £2 sin(y>)) 
H 3 y/W2 (sin(V>) (£3 csc(r?) - £4 cot(r?)) + £2 cos(ip)) 


n 1 ^ 


y/wiy/wi 

* 7 4 ^ 4 \/w 3 


y/my/W2y/W3 

y/Wly/W2 

1 

| 

l w\ 

! «4 , 


\Wi 

W 2 70 3 

y/W\W 2 W 3 ( 
A 

1 1 1\ 
„„2 + „„2 + ,„2 ) 


7 2 7 3 7 4 , 


where the following non-vanishing are used: 


e n = 


e, = 


y/W\W 2 W 3 ’ 


w\ csc(rj) cos(ip) 


y/W 2 W 3 

y/w^ cot(rj) sin (ip) 

y/W\W 3 


e l = - 


69 — 


y/wi sm(ip) 

y/w 2 w 3 

y/w^csc^r]) sin(^) 


y/W 1W3 


y/W 3 

y/W\W 2 


e 2 


Having the symbol of Z?, a direct calculation yields 


y/W2 cos (ip) 

y/W\W 3 

y/wi cot(p) cos(t/>) 
y/w 2 w 3 


cj(D 2 ){x 1 C) = p 2 (x,€) + + Po{x,£), 


where the expressions for the pk(x , £) are recorded below. The principal symbol of 
Z? 2 , which is homogeneous of order 2 in £, is given by 

p 2 (x,£) = ---cot 2 (?7) cos 2 (^) + ^3^ csc 2 (?7) cos 2 (^) 

WiW 2 W 3 

+^2U W 1 cot(77) sin(2'0) - t( 2 t; 3 wl csc(ry) sin(2^) 

—2^ 3 ^4 w\ cot(??) csc(p) cos 2 (ip) + £ 2 7e 2 sin 2 (7/;) 

+tp 2 w\ cot 2 (77) sin 2 (ip) — £2^4 w 2 cot(77) sin(2^) 

+^3^2 esc 2 (77) sin 2 (7/>) + P,2p,sw\ csc(?7) sin( 27 />) 

—2^3^4u;2 cot(?7) csc(p) sin 2 (7/;) + ^ 2,u, 2 cos 2 (ip) 

Hlwl+g)!, 
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where I is the 4 x 4 identity matrix. The component of cr{D 2 ) that is homogeneous 
of order 1 has a lengthy expression: 


Pi(x ,0 = 


(- 


it; 2 u>i cot(?7) cos 2 (ip) i^2'W2 cot(77) sin 2 (ip) 


W2W3 

3it;4W2 esc 2 (77) sin(2^>) 


wiw 3 

3*^4rui esc 2 (77) sin( 27 />) 


477:177:3 

7^47^2 cos(2?7) esc 2 (77) sin(2^) 
477:177:3 

7^377:2 cot(77) csc(r]) sin( 27 />) 


47772 7713 

7^477:1 cos( 27 y) esc 2 (77) sin( 2 ^>) 


+(^? + 


777 1 77 : 3 
^3 
2 t 7>2 


477:27773 

7^37771 cot(77) csc(r]) sin(2r/>) 

W2W3 


)I 


+( 


2 t 7 : 2 

7^47772 cot(77) sin(^) 


_ i^l) 7 2 7 3 
27773 n r 


i^3U>2 csc(r/) sin( 7 />) *£2^2 cos(t/>) 


+ 


277:4 

37772 csc(r/) sm(ip) 


27774 

7^47772 cot(77) sin(7/>) 

2777 § 

7^4 cot(77) sin(^>) 7^3 csc(rj) sin(f/>) 7^2 cos(^) 


27774 

it; 2 W 2 cos(t/>) 


277:3 


2^3 


2^,4 


+(- 


2w 2 

7^477:1 cot(?7) cos(t p ) 
2wl 


b i 


2w 2 2w 2 

i&wi csc(?7) cos(t/>) i&wi sin(^>) 


2 t 77 2 


2 t 77 § 


7^47771 COt(r]) COs(t/>) *£37774 CSc(r]) COs(t/>) sin( 7 /:) 


277:3 


2 t 7 :| 


2 t 77 § 


*£4 cot(77) COs(t/:) 


+(- 

+(- 


277:1 

7 ^ 4777 ^ 


*6 csc(?7) cos(V>) , 7^2 sin( 7 /») 3 4 

+ —^-b 7 


*£ 4«4 


277:1 

7^4777^ 


277:1 


1 _ 4 


■b 7' 


2 t 77 2 77’2 277:177:2 2W1W2W3' 

7^4 cot(77) COs(' 0 )t 77 , 1 7^3 csc( 77 ) COs(t/:) 777 i z £ 2 sin( 7 /)) 777 , 1 


+ 


+ 


277717772 7773 

7^4 cot(77) COs(V7)t7:2 
2777377:3 

7^4 cot (77) COs(i/))w ' 3 


2 w2Wo 


2W1W2W3 

it; 3 csc(?7) COs{'tp)w ' 2 
2777377:3 

*'£3 csc(?7) COs(lp)w 3 
27772 777 ? 


+ 


+ 


2W1W2W3 
ij2 sm(tp)w2 
2w 2 W3 

i £ 2 Sin(V:)777b ! 2 


277:277:? 


-bV 


+( 


*£4 cot(77) sin(7/>)7774 7^3 csc( 77 ) sin(^)) 777 i 7^2003(^)7774 


277747773 
7^4 cot(77) sin(7/:)TT72 
277:177727773 

*'£4 cot(77) sin(7/:)TT73 
2w\w 2 


277:477:3 

7^3 csc(rj) sm(%l))w2 
2W1W2W3 

7^3 csc(rj) sin(7/:)77:3 
277:477:3 


2w 2 U>3 
it;2 cos(i/j)w' 2 
277717772 7773 
*£2 COs(lp)w 3 
2wiw 2 


) 7 b 3 . 


7 



Finally we have the zero order part of cr(D 2 ): 

3 w[ 

8 wf 


Po(x,0 = (- 


w 


w\ 


Wo 


Wo 


Wo 


+ 


8W1W2 
3 w'o 


8w\w 2 


8 w 2 W 2 8 wfw 3 8 W 2 W 2 


8Wn 


(- 


+ 


W ' 3 + M) 7 3 4 + 

8 w 2 W 3 8 u >3 


W1W3 


^1^2 


■ + 


+ 


5w ' 2 


4w 2 W2W3 8w\w 2 W3 8wflU2W 2 16 wfw 2 W 3 


w' 2 w ' 3 


5 w! 2 


5 


8 w\w 2 w 3 

w 3 


I&W1W2W3 


w 2 


W1W3 


4W1W2W3 

W\ 


8w\W2 16 u >2 8W1W3 8W2W3 


16w 3 


)!■ 


4 w\W 2 W 3 

W2W3 


1GW1W2W3 

1171172. 


16 u;f 


3 . Calculation of the terms up to a 4 in the Spectral Action 

Calculation of the Seeley-de Witt coefficients associated with an elliptic positive 
differential operator on an m-dimensional compact manifold M can be achieved by 
using the Cauchy integral formula and parametric pseudodifferential calculus, see 
[ 26 ] for a detailed discussion. Let us review this method briefly for the operator D 2 , 
where D is the Dirac operator acting on a spin bundle on M. In order to derive a 
small time asymptotic expansion of the form 

OO 

Trace(e- tD2 ) ~ t m/2 o-2n(D 2 )t n (t ->■ 0+), 

n =0 

one can start with the Cauchy integral formula to write 

e ~ tD2 = — f e~ tx (D 2 - A)’ 1 d\ 

27 Tl J 1 

where the contour 7 in the complex plane goes around the non-negative real num¬ 
bers clockwise. Then, the idea is to approximate ( D 2 — A) -1 by pseudodifferential 
operators and to derive the above expansion by computing the trace of the corre¬ 
sponding approximation of the heat kernel. 

The symbol of D 2 is of the form p 2 {x^) +pi(a;, £) + po(a;,£) where each pi is 
homogeneous of order i in £. Since D 2 is an elliptic differential operator of order 2 , 
the inverse of D 2 — A is approximated by its parametrix R\ with 

OO 

<t(R x ) 

3=0 

where each 77 (x, £, A) is a parametric pseudodifferential symbol of order — 2 — j , in 
the sense that 

rj(x,t£,t 2 A) =t~ 2 ~ J r j (x,^,X). 

The equation 

can be solved recursively by comparing the homogeneous terms on the both sides 
after expanding the left hand side using the asymptotic composition rule for the 


a((D 2 - X)R x ) ~ ((p2(x,0 - A) + Pi(x,£) +p 0 (x,£)) o 


17=0 























symbols. Indeed, one finds that 


ro(x,€, A) = (P2(x,0 - A) \ 

and for any n > 1, the term r n (x,£, A) is found to be expressed in terms of 
ro(x, £, A),..., r n _i(a:, £, A) by the formula 

r n (x,£,\) = - ’S2—dg r j (x,£,\) D%p k (x,£)r 0 (x,t,\), 

where the summation is over all a G Z 4 0 , j G {0,1,..., n — 1 },k € {0,1,2}, such 
that |a| + j + 2 — k = n. 

The coefficients a 2n (D 2 ) in the small time asymptotic expansion are then com¬ 
putable by integrating the invariantly defined functions 

a 2n {x,D 2 ) = J J e~ x ti{r 2n (x,£ 1 A)) dAcP£ 

over the manifold against the volume form, which shows that these coefficients are 
local invariants of the geometry. We note that the odd coefficients vanish since for 
any odd j, the term rj(x,£, A) is an odd function of the variable £ € R m , whose 
integral over vanishes. 

Applying this method to Bianchi type-IX metric we compute the corresponding 
ao,a 2 ,a, 4 , which are recorded below (without writing the integral with respect t). 
The volume term is simply 

a 0 (D 2 ) = 4wiw 2 w 3 , 


and the scalar curvature term, after remarkable cancellations, is found to be 

K ) 2 K ) 2 


a 2 (D 2 ) = 


9 2 

wf w 2 2 

w l . 

xv\w\ 

2 2 

^1^3 , 

2 2 
U^Wg 

3 3 

3 + 

6 wl 

6 w 2 

Ow? 

K ) 2 

w[w 2 


xu 2 xjj' 3 

1 

6^1 

3w\w 2 

3w\W3 

3w 2 w 3 

3wi 


6 w 2 


6wf 


+ 


+ 


Although it seems lengthy, 
coefficient is expressed as: 


3w 2 3 w 3 

after an enormous amount of cancellations, the next 


«4 (D 2 ) 


W 1 W 2 W 1 W 3 

15w| “ 15w| 


wlwl 

15wf 


w\w 2 

15wf 


W\W 2 

15wf 


W?W 3 

15x4 


W 2 w 3 
15 w\ 


wiwf 

15x4 


W 2 Wg W\XJ0 2 W1W3 XJJ 2 XV 3 XJU 2 (w[) 2 XXJ3(xv[) 2 vo 3 (w' 2 ) 2 

1 " 1 - T ~T 7 - TT TT ^ - T" — 77 - T" "TT ^ 


15wf 15 w3 15w2 

Wl (w(,) 2 Wl (w() 2 

15wi 

W 2 (wg) 2 

15wiw| 15wiwf 15wfw 2 

2 (w () 2 | 2 (w 2 ) 2 

15w2w| 

15wyw3 

15w^W3 

I 5 W 1 W 2 W 3 I 5 W 1 W 2 W 3 

2 (W 3 ) 2 

w 2 (w {) 2 

w 3 (w }) 2 

Wl (w 2 ) 2 xv 3 (w 2 ) 2 

I 5 W 1 W 2 W 3 

18w^W3 

18w^W2 

18wfw3 18wiw| 

Wl (xv' 3 ) 2 

W 2 (wg) 2 

W 2 w 3 (w'l) 

2 W 1 W 3 (xu 2 ) 2 W 1 W 2 (wjj) 2 

18 W 2 W 3 

18wiw| 

18wf 

18wf 18wf 

31 (w {) 4 

31 (W 2) 4 

31 (w 3 

) 4 7w[xju 2 7w[w'3 7w 2 W3 


90wfw 2 W3 9 OXJU 1 XV 2 XV 3 90w\w 2 xv3 60wf 60 w 2 60wJ 
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w w 


1 UJ 2 


W 1 W 2 w' 2 w' 3 5 w3Vu[w' 2 


45 wfw 3 45 w 2 W 3 45u>iU>§ 36 wf 


5w3w[w 2 5w2W , 1 w' 3 


36wi 


36w^ 


5 w 2 w' 1 w 3 5wiw 2 w 3 5wiw 2 w 3 ^ 7 w 3 w' 1 w 2 ^ 7 w 2 'w' 1 w 3 ^ 7w\w' 2 w 3 


36 w% 
41K) : 


36 Wg 


90 wfw 2 90wfwg 99w 2 w 3 
41 K) 3 ^ 41u>j (w 3 ) 3 41 w' 2 {w' 3 ) : 


30>w 2 

41tui (w' 2 ) 3 

180wflU 2 W3 180wfw2'U’3 180wfu>2W3 180wiW 2 W 3 

23 K) 2 K) 2 23 K) 2 K) 2 23 K) 2 («,') 2 


180wfw 2 W3 
41(«4)V 3 
180uiiW2 w i 90 w 3 -u;2W3 

U>ju>g wjlTg W 2 W 3 

45u> 2 u>2 45w2W 3 47nu\w 2 

91w[w 2 (w 3 ) 2 W 2 W 1 W 3 W 1 
IZQwlwlwl + 24w| + 24w 3 
w'i 


90WiW2'W 3 

91 (w[) 2 w 2 w 3 
180w 3 u;|w| 


9 OW 1 W 2 W 3 


'l w 2 

.,// 


9lw[ [w' 2 ) 2 w ' 3 
\89w\w 2 w 3 


w 3 w 2 wiw 3 w 2 w 3 
24 iul + 24 w 3 2 + 24n; 3 


, 5 K) *vfj 

8 w\W 2 W 3 

71w 2 w 3 w 3 

vWs 

41 (w' 3 ) 2 w'i 

9 T 


W\W 2 

24w 3 

- + W3W1 W1W2 

36u> 2 W3 36w 2 W2 36w 2 W3 

" 5K)V | 5K)M 

8 wiW 2 W 3 

w" t 7 lw[w 3 W 3 
n "3 


W " W 2 w 3 W 2 Vj" W 3 I 

I 2 W 2 W 3 12w\W3 12w\W2 36w 2 W3 36 w- 
5 W 2 W 3 W 1 5 wiW 3 W 2 5 w\W 2 W 3 5 (w[) 2 w" 

72 wf 72 w 2 72w 3 + 8 wfw 2 W 3 

5 (wg) 2 w 3 71wiw' 2 w" 7\w'-yw 2 w 2 TluijwgW)" 71w[w 3 W3 

180wIw 2 W3 180 w\w\w3 180u> 3 u>2W§ 180 u> 2 u>2W§ 

, 71 | 41 (w 2 ) 2 w'[ | 41K) 2 < | 41( W ')M 

180wiw 2 w 3 360 w\w 2 W3 3Q9w\w2 , w 3 3601^1^2 ^3 

^- 3 ^- 2 . + 41K)X + 41 «)M + + 

360wiw|w| 3 Q 9 w\w 2 W 3 360w 3 w;2in§ 36w 2 w 2 w 3 36 w\w 2 w 2 

11 W' IW ' 2 W'3 _ «) 2 _ «) 2 _ K) 2 + W 3 w'j 
36w 2 w 2 w 2 6 w{W 2 W 3 6 wiw 2 W 3 6 wiW 2 W 3 3d>w\w 2 

W\W 3 W2W 3 w"w 2 W 2 W 3 vj"w 3 

36w2W 3 36wiw 3 15wfw 2 W3 15wi w 2 w 3 15wfw2iv 3 


w' lWl W 

6 wfw2W3 


w' 2 W2 {3) _ w' 3 W3 {3) _ w’ 2 W 1 < 3 > _ w' 3Wl ^ _ w[w 2 {3) _ w' 3 W 2 {3) 

Qw\W 2 W 3 6w\W2W 3 10wfw 2 W3 10wfw2W 3 10wfw 2 W3 10wiW 2 W 2 

w[w3^ w 2 w 3 ^ Wi (4) W 2 ^ W3^ 

10w 2 W2w| IOwi^Iw 2 30w 2 W2W3 30w\W 2 W3 30w\W2W 2 

Following the conjecture of Chamseddine and Connes for Robertson-Walker met¬ 
rics [S] , which was addressed in PH, the crucial observation to make at this stage is 
that all of the coefficients appearing in the above terms are rational numbers. This 
indicates that the rationality result holds for the Bianchi type-IX metric, which is 
proved in i}5] 


4. Heat Coefficients and the Wodzicki Residue 

In this section we introduce a method for computing the Seeley-de Witt coeffi¬ 
cients by making use of Wodzicki’s nonconnnutative residue @21148]. The advantage 
of this method is that it involves significantly less complexity in computations, thus, 
for instance, it illuminates the structure of the heat expansion of the Bianchi type- 
IX metric more elegantly. 
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Given a closed m-dimensional manifold M, the Wodzicki residue is the unique 
trace functional on the algebra of classical pseudodifferential operators acting on 
the smooth sections of a vector bundle over M (up to multiplication by a constant). 
The local symbol a of a classical pseudodifferential operator P a of order d £ Z has 
an asymptotic expansion of the form 


OO 


<7d-j(x,t) (£-too), 

3 = 0 


where each (Jd-j is positively homogeneous of order d—j in £. The noncommutative 
residue of the operator P CT is defined by 



where S*M = {(a;, £) £ T*M;\\^\\ g = 1} is the cosphere bundle of M and the 
integral is in fact the integral of the corresponding Wodzicki residue density over 
M, see [471 [481 [29] for more details. 

An alternative definition for Res, which is quite spectral, provides a link between 
the Seeley-de Witt coefficients and the noncommutative residue. That is, for any 
pseudodifferential operator P CT , the map that sends a complex number s with a large 
enough real part to Trace(P CT A _,s ), where A is a Laplacian, has a meromorphic 
extension to the complex plane with at most simple poles at its singularities. The 
noncommutative residue can be defined as the linear functional 


P(t i-t res s=0 Trace(P cr A s ), 


which turns out to be a trace functional. Thus, there is a constant c m such that 
for any classical P CT , we have 


Res(P CT ) = <J- m (x, £) d m d m x = c m (res s=0 Trace(P o -A s )) . 



im—1 


The constant c m can be computed easily as follows. The operator Res(A m / 2 ) is 
of order — m and its principal symbol is given by 


CTp ( A -m/2) = crp(A _1 ) m / 2 = crp(A) _m / 2 , 


which yields 


Res(A" m / 2 )= [ tr (a P (A)" m / 2 ) d m ~ 1 £,d m x. 


J S* M V J 

On the other hand, writing crp(A) = (^ i j /, we have 

res s= oTrace(A -m / 2 A -s ) = res s=m / 2 Trace(A _s ) 


• s—m 



r(m/2) 


li 








where for the second identity, we have used the formula @ proved below, for n — 0. 
Therefore, we have 


/ \ -m/2 

/(ZijMitj) <!'" 'Zd 


Cm = 2 m 7r m / 2 r(rn/2) 


= 2 


m+1 

7 r . 


f det~ 1 ^ 2 (g i ^) d m x 

Using the Mellin transform, A~ s = p^y / 0 °° e tx t s 1 dt, A > 0, one has 

1 Z 100 dt 

Trace(A _s ) = ——— / (Trace(e“* A ) — dim ker(A)) t s —. 

r(s) Jo t 

By breaking the interval of the integration in the latter to [0,1] and (1, oo), and by 
substituting the small time asymptotic expansion, 

N 

Trace(e" tA ) = t ~ m / 2 ^ a 2 n t n + 0{t ~ m / 2+N+1 ), 

n =0 


in the first part, one finds that 

res s=m/2 -„Trace(A~ s ) = Q2 »( A ) 

I (m/2 — n) 

for any non-negative integer n < m/2 — 1, see (26j . In particular we have 
(2) res s= iTrace(A” s ) = a m _ 2 (A). 

This observation yields the following assertion, which is used crucially in the sequel. 


Lemma 4.1. If A is a Laplacian acting on the smooth sections of a vector bundle 
over an m-dimensional manifold, then 

a m _ 2 (A) = — Res^ 1 ) = ^ Res(A~ 1 ). 

Cm 2 m+i 7T m 

Proof. It follows from the identity (0 and the fact that 

res s= iTrace(A~ s ) = res s= oTrace(A _1 A _s ) = —Res(A _1 ). 


Since we are mainly concerned with studying the spectral action for the Bianchi 
type-IX metric in this article, let us assume that D is the Dirac operator on a 
4-dinrensional manifold. By applying Lemma G3D to A = D 2 , we have 

a 2 (D 2 ) = —Res {D~ 2 ) = ( tr (<t_ 4 (D~ 2 )) d 3 £d 4 x, 

c 4 32?r 4 J S , M 

where cr_ 4 (IA~ 2 ) is the homogeneous component of order —4 in the expansion of the 
symbol of the parametrix of D 2 . In the following theorem, we show that the next 
coefficients a 2n (IA 2 ),n > 2, can similarly be expressed as noncommutative residues 
of Laplacians. 

Theorem 4.1. Let D be the Dirac operator on a f-dimensional manifold. For any 
non-negative even integer r, we have 

a ?+r(D 2 ) = 25 7r ! +r/2 Res(A- 1 ), 

where 

A = D 2 <g>l + l<g> Ar, 
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in which A-p- is the flat Laplacian on the r-dimensional torus T r = (R/Z) ? . 

Proof. It follows from Lemma 14.11 that 

a 2+r (A) = —Res(A _1 ), 

C4_)_ r 

where ct_ 4 _ r (A _1 ) is the homogeneous term of order —4 — r in the expansion of 
the symbol of the parametrix of A. 

Since the metric on T r is chosen to be flat, its volume term is evidently the 
only non-zero heat coefficient, which combined with the Kiinneth formula (cf. [,26j) 
implies that 

a 2 +r ((x,x') : A) = a 2 +r (x,D 2 )a 0 (x',T r ) 

= 2 ~ r 'K~ r ^ 2 a 2 +r {x,D 2 ). 

Therefore, 

or r / 2 

a 2 +r (D 2 ) = -Res(A _1 ) 

^4+r 

= 25^/2 ReS ( A " 1 )' 

□ 

A direct consequence of this theorem provides an efficient method for computing 
the Seeley-de Witt coefficients with significantly less complexities in the calcula¬ 
tions. It also yields an elegant proof of the rationality result for the Bianchi type-IX 
metric, which is presented in the following section. 

Corollary 4.1. Assuming the conditions and notations of Theorem EH we have 
a 2 +r (D 2 ) = 25 jT l +r/2 j tr(CT_ 4 -r(A- 1 )) d 3+r £ d 4 x. 

Proof. It follows from the fact that if <j(D 2 ) = p 2 (£,£) +pi(x, £) + po (£,£)< where 
each pi is homogeneous of order i in £, then cr(A) = p 2 (x, £) +Pi{x,£) +po(x,t;), 
with p 2 {x,f) = P 2 {x,%) + {f 2 + • • • + £, 2 +r )I- Therefore the homogeneous terms 
(A -1 ) °f or der —2 — j in the expansion cr(A _1 ) ~ Ujlo a -2-j (A' 1 ) are 
independent of the coordinates of T r . 

□ 

We confirm the validity of the coefficients ao,a 2 ,a 4 , calculated for the Bianchi 
type-IX metric in ij3] by noting that the method devised in the present section 
produces the same expressions. We stress that in practice the new method is 
significantly more convenient since the expression that leads to a Seeley-de Witt 
coefficient simplifies when one considers its restriction to the corresponding cosphere 
bundle in order to compute the noncommutative residue. 

The noncommutative residue was originally discovered in the 1-dinrensional case 
by Adler jTj and Manin [34] . Its coincidence with the Dixmier trace m on pseudo¬ 
differential operators of order —to on an m-dimensional closed manifold indicates 
its applicability for explicit and convenient computations. It is also worth men¬ 
tioning that a noncommutative residue developed for noncommutative tori |23| 122] 
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simplified a purely noncommutative heat kernel computation significantly and clar¬ 
ified in 13201 the reason for mysterious and remarkable cancellations that occur in 
this type of computations. 

5. Rationality of the Spectral Action for Bianchi Type-IX Metrics 

The Seeley-de Witt coefficients a 2n appearing in the expansion of the spectral 
action for the Bianchi type-IX metric are expressed in terms of several variable 
polynomials with rational coefficients evaluated on the cosmic evolution factors 
wi(t),w 2 (t),w 3 (t), and their derivatives of certain orders. This extends the state¬ 
ment conjectured in [9] and addressed in [21] for Robertson-Walker metrics to a 
homogeneous anisotropic cosmological model. 

In order to prove the rationality result for the Bianchi type-IX metric, similar to 
the treatment in [ 21 ], let us start with the crucial observation that the local forms 
a 2n (x, D 2 )d 3 x, where D is the Dirac operator of this geometry, are invariant over 
the spatial manifold § 3 . This can be seen from the defining formula m for the 
metric, in which the left invariance of the 1 -forms or, <r 2 , <r 3 , implies that the metric 
is invariant under any diffeomorphism arising from left multiplication by an element 
of SU{ 2 ). Since the action is transitive and left multiplication by an element of 
SU (2) is an isometry, any isometry-invariant function on § 3 is independent of the 
spatial coordinates. In particular, the restriction of the kernel of e~ tD to the 
diagonal and consequently the differential forms a 2 n (a;, D 2 ) dfx are invariant, and 
if we set a 2 n{x, D 2 ) drx = a 2n (x, D 2 ) dvol g , where dvol g is the volume form, then 
a 2n (: r, D 2 ) is independent of the spatial coordinates for any n. 

Furthermore, we can easily determine the general form of a 2 n (:r, D 2 ) by applying 
the method devised in [ 2 ] which is based on making use of the noncommutative 
residue, combined with the Kiinneth formula and restricting the computations to 
the cosphere bundle. In fact, writing cr(A _1 ) ~ a j ( x > £), where each Oj is 

homogeneous of order j, one finds recursively that 

v~ 2 (xffi) = p 2 (xffi)~ 1 , 

cr_ 2 -„(x,£) = -'^2^dga j (x,(;)D%p k (x,t;)(T- 2 (x,0 (n > 0), 

where the summation is over all multi-indices of non-negative integers a, —2 — n < 
j < — 2 , 0 < k < 2 , such that |a| — j — k = n. 

Thus, if we define C/i+i = Yv e u£H-ij then it can be shown by induction that 

(3) cr- 2 _ n (a;,£)ls*(MxT"-2) = a- 2 -n(x, C(0)l C eS"+i = (™i w 2 w?,)~i n P n {C), 

for any integer n > 2 , where P n ( C) is a polynomial in Ci>--',Cn+ 2 > with the 
coefficients being matrices whose entries are in the algebra generated by ratio¬ 
nal numbers, trigonometric functions of the spatial coordinates, and where 
i e {1, 2, 3}, p £ {0,1,..., n}. This fact leads to the following statement about the 
general form of the coefficients a 2 n (D 2 ). 

Theorem 5.1. For any non-negative integer n, the coefficient a 2n (Z) 2 ) in the ex¬ 
pansion of the spectral action for the Bianchi-type IX metric is of the form 

a 2 n[D ) = [ W1W2W3 ) Q 2 n \ Wl , W 2 , W 3 , W 1 , W 2 , W 3 ,..., W 1 , w' 2 ',u > 3 1 , 

where Q 2n is a polynomial with rational coefficients. 


14 


Proof. It follows from Corollary 14.11 that 

“ 2+ ' (i " d2) = 5v^X, < „ xT , ) tr( ‘ 7 - 1 - (A ' 1)) ‘ i3+ '« 

= l, tr dvol, d 3 +'(, 

where C/i+i = Yl v e u£w +1 so that the Jacobian of the coordinate transformation is 
just dvolg. This implies that 

a 2n {x,D 2 ) = 1 f tr(CT_ 2 - 2 n(A~ 1 ))d 2n+1 C, 

2 b Tt n+i J s2n+1 

which, as shown above, is independent of the spatial coordinates. Thus, we have 


02 n(D 2 ) 


a 2n {x, D 2 ) dvol 


g 


V0l-0 2 n{D 2 ) 
1 &T: 2 WiW 2 Wz CL2n{D 2 ) 


W 1W2W3 

2 7 r n + 1 


tr(<7_2-2ra(A -1 )) d 2n+1 (. 


The equation © allows us to write a- 2 - 2n (x, £(£)) = (w\w 2 w 3 ) 3 "P 2 ra (C), 
which yields 


a 2 n(D 2 ) = 


(wiw 2 w 3 ) 


1—3 n 


2 ir n+1 


tr(P 2n (C))(A- 1 )) d 2n+ \ 


{Wi W2W3) 1 3n Q 2 n (wi,W2,W 3 ,w[,W2,w' 3 , . . . , wf '" -) , 1 


(2 n) (2n)' 

^2 > w 3 


Note that tr (P2 n (0) is a polynomial in Ci; C 21 — > C 2 n+ 2 , with the coefficients in the 
algebra generated by the rational numbers, trigonometric functions of the spatial 
coordinates, and w where i £ {1, 2, 3}, p £ {0,1,, 2n}. 

The integral of a monomial m a {C) = c a C 1* 1 • ■ ■ C 2 n+ 2 2 over § 2n+1 is either 0, or 
can be written as 



m a (0 d 2n+1 C 


2c 0 n,r(^) 

T(n + 1 + ^) ’ 


if each aq is an even non-negative integer. Also, recall that T(^) = qir 2 for some 
q £ Q when n £ 2N + 1, and T©) £ Z when n S 2N. Therefore we have 


f m a (0d 2n+ \ = q^ =qn n +\ 

J S 2 "+1 

for some q £ Q if c a £ Q. Since 02 n(D 2 ) = ( wiw 2 w 3 ) 1 ~ 3n Q 2n has no spatial 
dependence, we conclude that 

Q 2 n = J tr(P 2n (C))(A" 1 )) d 2 ” +1 C 

belongs to the algebra generated by the and rational numbers. □ 
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6. Gravitational Instantons, Modular forms, and Rationality 


Among the Euclidean Bianchi type-IX models, an especially interesting class 
consists of the Bianchi IX gravitational instantons. A gravitational instanton is 
both self-dual (that is, the Weyl curvature tensor is self-dual) and an Einstein metric 
(the Ricci tensor is proportional to the metric). A remarkable feature of Bianchi 
IX gravitational instantons with SU (2) symmetry is that they can be completely 
classified in terms of solutions to Painleve VI integrable systems, [271 HU[46]. The 
latter are a 4-parameter family of singular ordinary differential equations of the 
form 


d 2 X 
dt 2 


1 


+ 


1 


+ 


1 


X X-l x-t 
X(X-1 )(X-t) 


t 2 (t- 1) 2 


dX 

dt 

t 


t - 1 X-t 
, **(*-1) 


a + p x 2 +1 (x- 1 ) 2 ww (x-ty 


dX 

dt 


The self-dual equation for the 517(2) Bianchi IX metrics is written in [33] as 
an ordinary differential equation in the Wi and in additional functions a,, i = 
1,2,3 that arise as the components of the connection 1-form in a basis of anti- 
self-dual 2-forms, see [461 . In terms of the conformally invariant variable x = 
(a 2 — cii)(a 2 — CI 3 ) 1 the self-dual equations for the Riemannian Bianchi IX metric 
can be rephrased as a system of equations 

Wi = Clix'(x(l — x )) _1//2 , 


n; = -■ 


^2^3 


O' — — 
, 5 ‘2 — 




it, = - 


fiiri2 


x(l — x) Z X 3 1 — X 

These in turn can then be reduced to a case of the Poincare VI equation with 
parameters 

\ a > PilA) (g) gi g) g)> 

see [331 US]. In [ 3 ], the solutions to this equation are given explicitly in terms of a 
parameterization involving theta functions and theta characteristics 

it\p,q]{z,in) ■= E exp (—7T (m + p) 2 p + 2 ni(m + p)(z + q )) . 

mGZ 

Namely, with the notation d[p,q\ := g]( 0 ,i/z), and 

i? 2 := 0 [ 1 / 2 ,O], 0 3 := 0 [O,O], — ^[ 0 ,1/2], 

one finds ai = 2 log 1 and 

* , , + _ i , , ^% + jg + g] 

Wl 2^ 3 4 e lrip, d[p, q\ ’ W2 2^ 2 4 e ni H\p,q] 

1 0 Q 3 ;% +id 
w3= %, ? i - 

The asymptotics of these solutions were analyzed in 135 ) . where it is shown that, 
for large /z, they approximate Eguchi-Hanson type gravitational instantons with 

W2 = W 3 ^ W \, un¬ 
it is clear that, for the Bianchi IX gravitational instantons, using the parameter¬ 
ization of [ 3 ], the Seeley-de Witt coefficients 02™ of the spectral action are rational 

16 
















functions, with Q-coefRcients, in the $2, $3, $4, i?[p, 5], < 9 g $[p, <7] and e l7rp and deriva¬ 
tives, hence they belong to the field generated, over Q, by these functions. We will 
return in a second part of this work [ 19 ] to discuss in detail the arithmetic properties 
of the spectral action for Bianchi IX gravitational instantons. 

In this perspective, one can view the rationality question about the spectral ac¬ 
tion in a similar light to analogous questions that occur whenever arithmetic and 
number theoretic structures arise in theoretical physics. For example, when Feyn¬ 
man integrals are interpreted as periods (see [ 35 ] for an overview of that setting), 
the fact that the relevant amplitude forms and domains of integration are algebraic 
over Q (or Z) has direct implications on the class of numbers that arise as periods. 
Another such instance of arithmetic structures in physics, where rational coeffi¬ 
cients play an important role, is in the zero temperature KMS states of quantum 
statistical mechanical systems: in the case constructed in m (see also Chapter 
3 of jT 5 j) for instance, the construction of an arithmetic algebra of observables, 
defined over Q, is linked to modular functions and makes it possible to have KMS 
states with values in the modular field. The relation between the spectral action of 
Bianchi IX gravitational instantons and modular forms will be discussed in m- 


7. Conclusions 

We have shown that the Seeley-de Witt coefficients a 2n (.D 2 ) associated with the 
Dirac operator D of the Bianchi type-IX metric, which appear in the expansion of 
the spectral action [ 5 ] , are expressed by polynomials with rational coefficients eval¬ 
uated on the cosmic evolution factors wi(t),W2(t),W3(t), and their derivatives of 
certain orders. It is quite interesting that although this metric provides a homoge¬ 
neous anisotropic cosmological model, after remarkable cancellations, only rational 
coefficients appear in the final expression for each a,2 n (D 2 ). Such a rationality re¬ 
sult was first conjectured in [ 9 ] for Robertson-Walker metrics, which was addressed 
in [2T| . 

Our proof of the rationality statement for the Bianchi type-IX model, similar to 
the argument given in [21] . begins with the crucial observation that the kernel of 

1 T-\ 2 

e is restricted to have no spatial dependence on the diagonal. We then take a 
novel approach to proceed the argument. That is, we have devised a general method 
that expresses the Seeley-de Witt coefficients of a geometry as noncommutative 
residues of operators. This is an efficient method that allows explicit calculations 
with significantly less complexities, compared to the method of using parametric 
pseudodifferential calculus [25 . More importantly, it leads to an elegant proof of the 
rationality result for the Bianchi type-IX metric. To be more explicit, the Wodzicki 
residue 03S5I involves an integration over the cosphere bundle of a manifold, and 
the expression for computing a 2 n{D 2 ) simplifies to our favor when restricted to the 
cosphere bundle, in the view of our method. 
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